Interacting vortex pairs in inviscid and viscous planar flows 



Thierry Gallay 

Universite de Grenoble I 
Institut Fourier, UMR CNRS 5582 
B.P. 74 

F-38402 Saint-Martin-d'Heres, France 
Thierry . Gallay@uj f -grenoble . f r 



o 

(N 

O 
D 

Q I December 10, 2012 

^ —| " Abstract 

cm: 

<^ The aim of this contribution is to make a connection between two recent results concerning 

the dynamics of vortices in incompressible planar flows. The first one is an asymptotic 
expansion, in the vanishing viscosity limit, of the solution of the two-dimensional Navier- 
Stokes equation with point vortices as initial data. In such a situation, it is known [5] that 
the solution behaves to leading order like a linear superposition of Oseen vortices whose 
centers evolve according to the point vortex system, but higher order corrections can also 
be computed which describe the deformation of the vortex cores due to mutual interactions. 
1 The second result is the construction by D. Smets and J. van Schaftingen of "desingularized" 

solutions of the two-dimensional Euler equation [22]. These solutions are stationary in a 
uniformly rotating or translating frame, and converge either to a single vortex or to a vortex 
pair as the size parameter e goes to zero. We consider here the particular case of a pair 
of identical vortices, and we show that the solution of the weakly viscous Navier-Stokes 
equation is accurately described at time t by an approximate steady state of the rotating 
Euler equation which is a desingularized solution in the sense of [22] with Gaussian profile 
and size e = yfut. 
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1 Introduction 



Numerical simulations of freely decaying turbulence show that vortex interactions play a crucial 
role in the dynamics of two-dimensional viscous flows [12, 13]. In particular, vortex mergers 
are responsible for the appearance of larger and larger structures in the flow, a process which 
is directly related to the celebrated "inverse energy cascade" [3]. Although nonperturbative 
interactions such as vortex mergers are extremely complex and desperately hard to analyze 
from a mathematical point of view [11, 20], rigorous results can be obtained in the perturbative 
regime where the distances between the vortex centers are large compared to the typical size of 
the vortex cores. 

As an example of this situation, consider the case where the initial flow is a superposition of 
./V point vortices. This means that the initial vorticity loq satisfies 

N 

= ^aidi- - Xi) , (1.1) 

i=i 

where x\,...,xjq S M 2 are the initial positions and ai,...,ajv £ K the circulations of the 
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Figure 1: The motion of two point vortices whith circulations a\ > a 2 > (left) and ol\ +a 2 = (right). 

vortices. Let u(x,t) be the solution of the two-dimensional vorticity equation 

+ u-Vuj = uAu , xeR 2 , t>0, (1.2) 

at 

with initial data ojq, where u(x,t) is the velocity field defined by the Biot-Savart law 

"(M) = ^- / ^%w(j/,t)dy , xGM 2 , t>0. (1.3) 
2vr 7 R 2 |x - y\ 2 

Solutions of (1.2), (1.3) with singular initial data of the form (1.1) were first constructed by 
Benfatto, Esposito, and Pulvirenti [2]. More generally, if ojq G is any finite measure, 

Giga, Miyakawa, and Osada [8] have shown that the vorticity equation (1.2) has a global solution 
with initial data ujq, which moreover is unique if the total variation norm of atomic part of luq 
is small compared to the kinematic viscosity v. This last restriction has been removed recently 
by I. Gallagher and the author [4], so we know that (1.2) has a unique global solution to G 
C°((0, oo), L 1 (M 2 )nL 00 (lR 2 )) with initial data (1.1), no matter how small the viscosity coefficient 
is. This solution is uniformly bounded in L 1 (R 2 ), and the total circulation f R2 u(x,t) dx is 
conserved. 

In the vanishing viscosity limit, the motion of point vortices in the plane is described by 
a system of ordinary differential equations introduced by Helmholtz [9] and Kirchhoff [10]. If 
zi(t), . . . , Z]\r(t) G M 2 denote the positions of the vortices, the system reads 

d z(t) - 1 y- a (*(*)- sift))- 1 i - 1 N (14) 

and the initial conditions Zj(0) = X{ for i = 1, ... ,7V are determined by (1.1). A lot is known 
about the dynamics of the point vortex system (1.4), see e.g. [21] for a recent monograph devoted 
to this problem. Most remarkably, (1.4) is a Hamiltonian system with TV degrees of freedhom, 
which always possesses three independent involutive first integrals. In particular, system (1.4) is 
integrable if TV < 3, whatever the vortex circulations a±, . . . , may be. In the simple situation 
where N = 2, both vortices rotate with constant angular speeed around the common vorticity 
center, see Fig. 1 (left). In the exceptional case where ot\ + «2 = 0, there is no center of vorticity 
and the vortices move with constant speed along parallel straight lines, see Fig. 1 (right). 
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It should be remarked that system (1.4) is not always globally well-posed: if N > 3 and if the 
circulations a±, . . . , ajv axe not all of the same sign, vortex collisions may occur in finite time for 
exceptional initial configurations [19, 21]. In what follows, we always assume that system (1.4) 
is well-posed on some time interval [0, T], and we denote by d the minimal distance between any 
two vortices on this interval: 



d = min min \zi(t) — Zj(t)\ > . (1.5) 
te[o,T] i+j 3 

We also introduce the turnover time 

d 2 

Tq = - — - , where \a\ = \a\\ + • • • + \cxn\ , (1.6) 
\a\ 

which is a natural time scale for the inviscid dynamics described by (1.4). For instance, for a 
pair of vortices with circulations of the same sign, one can check that the rotation period of each 
vortex around the center is 47r 2 T . 

When the viscosity v is nonzero, the point vortices in the initial data (1.1) are smoothed out 
by diffusion, and the solution uj(x,t) of (1.2) is no longer described by the point vortex system. 
In the particular case where N = 1, the unique solution of (1.2) with initial data ujq = ado is 
the Lamb- O seen vortex: 



' vt ' \Jvt \\Jvt' 

where the vorticity profile G and the velocity profile v G have the following explicit expressions: 

« = r H!|,,, l ^-iw( 1 - e ~ mi )- (L8) 

As was shown by C.E. Wayne and the author, Oseen vortices describe the long-time asymptotics 
of all solutions of the two-dimensional Navier-Stokes equation for which the vorticity distribution 
is integrable, and are also the only self-similar solutions of this equation with integrable vorticity 
profile [7]. 

When iV > 2, the solution of (1.2) with initial data (1.1) is not explicit, but in some 
parameter regimes it can be approximated by a linear superposition of Oseen vortices whose 
centers evolve according to the point vortex system (1.4). More precisely, we have the following 
result: 



Theorem 1.1 [5] Given pairwise distinct initial positions x±, . . . ,xn G 1R 2 and nonzero circu- 
lations a±, . . . , ajv £ fix T > such that the point vortex system (1.4) is well-posed on the 
time interval [0, T]. Let d > be the minimal distance (1.5) and Tq the turnover time (1.6). 
Then the (unique) solution of the two-dimensional vorticity equation (1.2) with initial data (1.1) 
satisfies 

±[ u;(x,t)-J2^G(^M) dx<K^ 2 , te(0,T], (1.9) 
\ a \ Jk 2 t^i 

where z{t) = {zi(t), . . . , zjy(t)} is the solution of (1.4) and K is a (dimensionless) constant 
depending only on the ratio T/Tq. 

Theorem 1.1 gives nontrivial information on the solution of (1.2) in the weak interaction 
regime, where the size 0(y/ui) of the vortex cores is much smaller than the distance d between 
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the centers. If the initial data are fixed, a convenient way to achieve this is to assume that the 
viscosity v is small compared to the total circulation |a|, and that the observation time T is 
smaller than or comparable to the turnover time To. Indeed, using definition (1.6), we see that 

vt ut 
d? \a\ T 

In particular, in the vanishing viscosity limit, we obtain the following corollary. 

Corollary 1.2 [5] Under the assumptions of Theorem 1.1, the solution oj u (x,t) of the viscous 
vorticity equation (1.2) with initial data (1.1) satisfies 

N 

w v {-,t) » V <*<*(•-*(*)) , for all t G [0, T] , (1.10) 

v— H) ^ — * 

i=l 

where z(t) = {z±(t), . . . , 2jv(i)} is the solution of (1.4). 

In other words, the solution of the vorticity equation (1.2) with initial data (1.1) converges 
weakly, in the vanishing viscosity limit, to a superposition of point vortices which evolve ac- 
cording to the point vortex dynamics (1.4). In particular, Corollary 1.2 provides a natural and 
rigorous derivation of the point vortex system itself from the Navier-Stokes equation. As is 
well-known, system (1.4) can also be rigorously derived from Euler's equation [15, 18], but in 
the latter approach it is necessary to regularize the initial vorticity because we do not know how 
to solve Euler's equation with singular data such as (1.1). In this respect, it is important to 
keep in mind that the right-hand side of (1.10) is not a weak solution of the inviscid vorticity 
equation dfOJ + u ■ Vw = 0. For completeness, we also mention that the vanishing viscosity 
limit for solutions of (1.2) with concentrated vorticity has been studied by Marchioro [16], who 
obtained the analog of Corollary 1.2 in that context. 

The conclusion of Theorem 1.1 can be interpreted in the following simple and somewhat 
naive way: If we solve the two-dimensional vorticity equation (1.2) with point vortices as initial 
data, the diffusion term uAoj in (1.2) smooths out the point vortices into Oseen vortices, and 
the advection term u ■ Vw translates the vortex centers according to the point vortex dynamics 
(1.4). While correct, this interpretation ignores the important fact that the advection of a 
smooth vortex by the inhomogeneous velocity field created by the other vortices results non 
only in a translation of the vortex center, but also in a deformation of the vortex core. In our 
case, this deformation is of order 0(vt/d 2 ) in L 1 norm, and therefore does not appear in (1.9) 
because it is included in the error term. However, as is explained in [5], such a small deformation 
of the vortex profile creates a self-interaction effect of size 0(1), which basically counterbalances 
the influence of the external velocity field, except for a rigid translation. 

As a matter of fact, self- interactions play a crucial role in the proof of Theorem 1.1, and 
we even believe that is it not possible to establish (1.9) without computing a higher order 
approximation of the solution. A systematic asymptotic expansion is carried out in [5], but if we 
only keep the leading order nonradial corrections to the Oseen vortices we obtain the following 
approximate solution of (1.2): 

where 

*«. o - «*) e 5 j^^^j^w - + °(r) ■ ^ 
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Figure 2: The level lines of the vorticity distribution for a pair of vortices with equal (left) or opposite 
(right) circulations. 



Here a : (0, oo) — > R is a smooth, positive function satisfying a(r) ~ C\r 2 as r — > and 
a(r) ~ (72r 4 e _r / 4 as r — > oo for some C\,C2 > 0. In the regime where the viscosity v is much 
smaller than the circulations Oj of the vortices, we can use formulas (1.11), (1.12) to compute, 
to leading order in our expansion parameter ut/d 2 , the deformation of the vortex cores due to 
mutual interactions. For any fixed t £ (0, T], this first order correction depends only on the 
relative positions of the vortex centers, which are determined by (1.4). In the particular case of 
a pair of vortices with equal or opposite circulations, the level lines of the vorticity distribution 
cj ap p(x,£) are represented in Fig. 2. 

The discussion above shows that the dynamics of weakly interacting viscous vortices is es- 
sentially driven by two different mechanisms: diffusion, which is responsible for the continuous 
growth of the vortex cores, and advection, which creates the motion of the vortex centers and the 
deformation of the vortex profiles. The latter effect persists in the vanishing viscosity limit, and 
it is therefore reasonable to expect that, if we can find solutions of Euler's equation describing 
widely separated Gaussian vortices, these inviscid solutions will provide an accurate approxima- 
tion of the viscous N- vortex solution considered in Theorem 1.1, if f is sufficiently small. The 
aim of this contribution is to explore this idea in the particular case of a single vortex pair. In 
this simple situation, we only need to consider solutions of Euler's equation that are stationary 
in a uniformly rotating or translating frame. 

Assume thus that N = 2 and, for definiteness, that both circulations 0:1,0:2 are positive. 
Given n, r2 > such that oiri = 0:2^, let d = r\ + r2 and Q, = (01 + 02) /{2nd 2 ). We consider 
the inviscid vorticity equation in a rotating frame with angular speed : 

d t oj + (it- Qx 1 ) ■ = . (1.13) 

Formally, the vorticity distribution 



ujq = oi<5(- — x\) + 02<J(- — X2) , where x\ = 




(1.14) 



is a stationary solution of (1.13). In the laboratory frame, this corresponds to a time periodic 
solution of (1.4) where both vortices rotate around the origin with angular velocity J). Now, 
let w* G 5(R 2 ) be a nonnegative vorticity profile which is radially symmetric, decreasing along 
rays, and normalized in the sense that J R2 w*(x) dx = 1. Given e > 0, we look for a stationary 
solution of (1.13) of the form 
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where Xi >t — > X{ and Wi t€ — > w* as e — > 0. Existence of such " desingularized" solutions of 
Euler's equation has been investigated in a recent work by D. Smets and J. Van Schaftingen 
[22]. In fact, the authors of [22] do not consider rotating vortex pairs of the form (1.15), but 
they treat a variety of other interesting cases, including a single stationary vortex in a bounded 
or unbounded domain, a rotating vortex in a disk, and a translating vortex pair in the plane. 
They use the "stream function method" , which consists in constructing (by variational methods) 
a nontrivial solution to an elliptic equation of the form — Aip = f e (ip + ^J7|x| 2 ), where f e is a 
power- like nonlinearity which depends on e in an appropriate way. The vorticity to = —Aip is 
then a stationary solution of (1.13). For technical reasons, the desingularized vorticity profiles 
obtained in [22] are always compactly supported, but it is reasonable to expect that similar 
results can be obtained with Gaussian profiles too. We hope to clarify this issue and to extend 
the results of [22] to pairs of vortices of the same sign in a future work. 

In Section 2 below, we study in detail the case of two identical vortices (a\ = 02). For a 
large class of radially symmetric profiles w*, we prove the existence of approximate stationary 
solutions of (1.13) of the form (1.15). In other words, we construct an asymptotic expansion 
in powers of e of the vorticity distribution (1.15) as a steady state of (1.13). Under natural 
symmetry assumptions, we show that this expansion can be performed to arbitrarily high order. 
Then, in Section 3, we prove that the solution u: u (x,t) of the rotating viscous vorticity equation 
dtU) + (u — Qx- 1 -) ■ Vw = uAlo with initial data (1.14) is very close to the inviscid stationary 
solution lo 6 with asymptotic profile w* = G, if e = \fvt is sufficiently small. This means 
that, when the viscosity is small, the solution of (1.2) with initial data (1.14) slowly travels 
through a family of uniformly rotating solutions of Euler's equation, whose vorticity profiles are 
approximately Gaussian and evolve diffusively. We expect a similar picture to be relevant in the 
general situation considered in Theorem 1.1, although the corresponding inviscid solutions may 
be more difficult to identify in that case. 

Acknowledgements. I am grateful to Felix Otto for suggesting the point of view adopted in 
this paper, and to Didier Smets for many stimulating discussions. 

2 Approximate steady states of Euler's equation 

The aim of this section is to construct an asymptotic expansion for a family of stationary 
solutions of the inviscid vorticity equation, which correspond to weakly interacting vortex pairs. 
For simplicity, we only consider the particular case where both vortices have the same circulation 
a > 0. As is explained in the introduction, weak interaction means that the distance d between 
the vortex centers is large compared to the size of the vortex cores. If the vorticity distribution 
is given by (1.15), this condition is satisfied if e > is sufficiently small. Here, we find it more 
convenient to fix the size of the vortex cores, and to assume that the distance d between the 
centers is large. This alternative point of view is of course equivalent to the previous one, up 
to a rescaling. Note, however, that the rotation speed $7 will now depend on d and behave like 
a/ (nd 2 ) as d — > 00. 

From now on, we fix a > 0, d 3> 1, and we look for a stationary solution oj of (1.13) describing 
a pair of identical vortices with circulation a. We make the ansatz 

uj(x) = aw{x — Xd) + aw (— x — x^j , u(x) = av(x — Xd) — av(—x — Xd) , (2.1) 

where Xd = (d/2,0), w is a localized vorticity profile to be determined, and v = K[w] is the 
velocity field obtained from w via the Biot-Savart law (1.3). We assume that w belongs to the 
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Schwartz class 5(M 2 ), is nonnegative, and satisfies the normalization condition 

/ w{x)dx = 1 . (2.2) 
Jr 2 

We also impose the following symmetry 

w(x 1 ,-x 2 ) = w(x 1 ,x 2 ) , v 1 (xi,-x 2 ) = -vi(xi,x 2 ) , V 2 (X!,-X 2 ) = v 2 (xi,x 2 ) , (2.3) 

which implies that uj{—X\,x 2 ) = oj(x\,—x 2 ) = uj{xi,x 2 ) for all x = {x\,x 2 ) G M 2 . Finally, 
without loss of generality, we assume that 

/ x 1 w(x)dx = . (2.4) 

JR 2 

This means that the vorticity distribution uj defined in (2.1) is indeed a superposition of two 
localized vortices centered at the points ±x d - 

The distribution uj will be a stationary solution of the rotating vorticity equation (1.13) if 
the profile w satisfies 

(av(x — Xd) — av(—x — Xd) — Ox -1 ^ • Vw(x — Xd) = , x£l 2 . 

Replacing x by x + Xd and denoting Q = Q/a, we obtain the equivalent equation 

(v(x) - v(-x - 2x d ) - f2(x + Xd) 1 ^ ■ Vw(x) = , xGR 2 . (2.5) 

Note that the rotating term Q(x + Xd) ± behaves like the velocity field v(x) when x 2 is changed 
into — x 2 ; this implies that the symmetry (2.3) is indeed compatible with Eq. (2.5). To determine 
the rotation speed we multiply both members of (2.5) by x 2 and we integrate by parts over 
M. 2 . We easily obtain 

Q / (xi + d/2)w(x) dx = / iv 2 (x)—v 2 (—x — 2xd))w(x)dx. 
Jr 2 Jr 2 ^ ' 

This identity can be simplified if we use (2.2), (2.4), and the fact that J R2 v 2 wdx = 0. We thus 
arrive at the following relation, 

n = Q[ w ] ■= -- I v 2 {-x - 2x d )w(x) dx , (2.6) 
d Jr 2 

which determines as a function of w. 

As we shall see, since the vorticity profile w is normalized by (2.2), the corresponding velocity 
field v = K[w] satisfies 

-v 2 (-x - 2x d ) ~ , as d^oo, 

for any fixed x £ M 2 . In view of (2.6), this means that $1 ~ (vrti 2 )" 1 as d — > oo. In particular, if 
we take formally the limit d — > oo in (2.5), we obtain the limiting equation v ■ Vw = 0, which 
simply means that w (or v) is a stationary solution of Euler's equation. In what follows, we 
assume that the limiting profile w* is radially symmetric and stable in the sense of Arnold [1, 19]. 
Roughly speaking, this means that w*(x) is a strictly decreasing function of |x|. Given any such 
profile, we shall construct perturbatively a family of approximate solutions of (2.5), indexed 
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by the parameter d > 1, which converge to w* as d — > oo. Under natural assumptions, these 
approximate solutions are uniquely determined by the asymptotic profile w*. 

An important question that we leave open here is whether we can actually construct exact 
solutions of (2.5) which converge to w* as d — > oo, in which case our approximate solutions could 
be recovered by truncating the asymptotic expansion of the exact solutions. As was mentioned 
in the introduction, it should be possible to prove the existence of such solutions, at least for 
a particular class of compactly supported profiles by adapting the variational techniques of 
D. Smets and J. van Schaftingen [22]. It might also be possible to construct exact solutions of 
(2.5) for more general profiles using a fixed point argument of Nash-Moser type. We hope to 
clarify these issues in a future work. 



2.1 Asymptotic profile and functional setting 

Let w* G iS(M 2 ) be a radially symmetric, nonnegative function satisfying the normalization 
condition (2.2), and let = AT[it>*] be the velocity field obtained from w* via the Biot-Savart 
law (1.3). We can thus write 

Wm (x) = -q(\x\ 2 ) , v*(x) = ^-^7oQ{\x\ 2 ) , x£K 2 , (2.7) 
IT 2tt 

where q : [0, oo) — > R + is a smooth, rapidly decreasing function and Q(r) = J^q(s)ds. Note 
that Q(r) — > 1 as r — > oo. We assume that w* satisfies the following strong stability conditions : 

2 I ( \ 

q(r) < for all r > , and sup < 1 . (2.8) 

r>0 Q{r) 

The first condition in (2.8) implies of course that q(r) > for all r > 0, so that Q(r) > for all 
r > 0. In particular, compactly supported asymptotic profiles w* are excluded. This condition 
also implies that w* is a stable solution of the two-dimensional inviscid vorticity equation, with 
respect to perturbations in L 1 n L°° [17]. The second assumption in (2.8) is more technical in 
nature, and can probably be relaxed. It is satisfied, for instance, if q(r) = 7e~ 7r for some 7 > 0. 
Note that we always have 

-r 2 q'(r) 1 
sup — > - . 

r>o Q(r) 4 

Indeed, if we assume on the contrary that Q(r) + 4r 2 q'(r) > for all r > 0, then the function 
h(r) = Q"(r) + Q(r)/(4r 2 ) is nonnegative and satisfies h(r) ~ q(0)/(4r) as r — > 0, and h(r) ~ 
l/(4r 2 ) as r — > 00. Thus yjrh G L 1 ((0, 00)), but if we integrate by parts we obtain 

^ Vr~h(r) dr = ^(q"{t) + ^Q{r)) dr = , 

which yields a contradiction. Finally, in addition to (2.8), we also assume that q 2 /q' decays 
rapidly at infinity : 

k ( \2 

sup -. q ,Yl < 00 , for all k G N . (2.9) 
r>o \q'(r)\ 

As was already observed, the asymptotic profile w* is already an approximate solution of 
(2.5), (2.6) in the sense that, if we substitute {w* , ) for (w,v) in (2.5), the left-hand side 
converges to zero as d — > 00. Our goal is to construct here more accurate approximations, which 
take into account the interaction of the vortices. We look for solutions of the form 

w = w* + uj , v = v* + u , (2.10) 
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where u = K[ui] is the velocity field obtained from uj via the Biot-Savart law (1.3). The symmetry 
(2.3) implies that 

uj(xi,-x 2 ) = uj(xi,x 2 ) , ui(xi,-x 2 ) = -ui(xi,x 2 ) , u 2 (xi,-x 2 ) = u 2 (xi,x 2 ) , (2.11) 
and in agreement with (2.2), (2.4) we impose 

/ uj(x)dx = / x 1 uj(x)dx = . (2.12) 
Jr 2 Jm. 2 

Finally, we assume without loss of generality that uj has no radially symmetric component, 
namely 

i-2-k 

I uj{r cos 6,r sin 9) d6 = for all r>0. (2.13) 
Jo 

We can always realize (2.13) by including, if necessary, the radially symmetric part of uj into the 
asymptotic profile u>*. In this respect, it is important to note that both conditions in (2.8) are 
open. 

Inserting (2.10) into (2.5), we obtain for uj the following equation 

Auj + u ■ Vuj + R d [uj] = , (2.14) 
where A is the linearized operator defined by 

Auj = • Vuj + u • Vw, , (2-15) 
and Rd[uj] is a remainder term which depends on the distance d between the vortex centers 

R d [uj](x) = (v*(x + 2x d ) - u(-x - 2x d ) - &{w* + cu](x + Xd) 1 -^ ■ V(u>*(x) + uj(x)) . (2.16) 

In (2.14)-(2.16), it is understood that u = K[uj] is the velocity field associated to uj. 
We look for solutions uj of (2.14) in the Hilbert space 

I = (weL 2 (K 2 ) f \uj(x)\ 2 p(\x\ 2 )dx < oo) , (2.17) 

equipped with the scalar product 

(uji,uj 2 ) = / uji(x)uj 2 (x) p(\x\ 2 ) dx , uji,uj 2 e X , 
Jr 2 

and with the associated norm ||w|| = {uj,uj) 1 / 2 . Here the weight p : [0, oo) — > R + is defined by 

P(r) = ~TF~\ ' r > • ( 2 - 18 ) 
q'{r) 

The reason for this particular choice is that the linear operator A has nice properties in the 
space X, see Section 2.2 below. In view of (2.9), the asymptotic profile u>* and all its moments 
belongs to X : for any k G N, we have 



\x\ 2k wJ\ 2 



If 1 f°° 

— ^ / \x\ Ah q(\x\ 2 ) 2 p(\x\ 2 ) dx = - r 2k q(r) 2 p(r)dr < oo . 
Jr 2 Jo 
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It is easy to verify that the operator A commutes with the rotations about the origin in R 2 , 
see Lemma 2.2 below. It is thus natural to use polar coordinates (r, 9) in the plane, and to 
decompose our space X as a direct sum 

oo oo 

x = e x n = e p n x , (2.19) 

n=0 n=0 

where P n is the orthogonal projection in X defined by the formula 

2 x r^is 

(P n oj)(r cos9,rsin9) = — / uj(r cos 9', r sin#') cos(n(9 — 9')) d9' , nGN. 

In particular, Xq = PqX is the subspace of all radially symmetric functions, and for n > 1 
the subspace X n = P n X contains functions of the form uj(r cos 9, r sin 9) = ai(r)cos(n#) + 
a2(r) sin(n#). With this notation, condition (2.13) means that Pqoj = 0. 

2.2 The linearized operator and its right-inverse 

We now discuss the main properties of the linearized operator A defined in (2.15). In the 
particular case where w* is the profile G of Oseen's vortex (1.8), the operator A was studied in 
detail in [7, 14], and we shall obtain here analogous results in a more general situation. 

From (2.15) we know that A = Ai + A2, where Aioj = v* ■ Vw and A20J = K[oj] ■ Vu>*. As 
is easily verified, A2 is compact in X, while Ai is unbounded. The maximal domain of A is 
therefore 

D(A) = D(Ai) = {u G X I ■ Voj G X} . (2.20) 
The most remarkable property of this operator is that it is skew- symmetric in X. 

Lemma 2.1 For allcoi,oj2 G D(X), we have (Awi,^) + (wi,Acj2) = 0. 

Proof. We shall prove in fact that both operators Ai, A2 are skew-symmetric. First, since the 
weight p(|x| 2 ) is radially symmetric, we have 



{v* 



Vwi,w 2 ) + (wi,v* • V0J2) = / p(\x\ 2 ) ■ V (UJ1UJ2) dx = 

Jr 2 



because the velocity field p(\x\ 2 )v*(x) is divergence-free. Next, since Vw»(x) = (2/ir)xq' (\x\ 2 ), 
we have 

(Ui ■ V«J»,W 2 ) + (Wl,«2 • Vtt7*) = -— / ((x ■ Ui)0J 2 + (x ■ U2)0J\) dx = , 

see e.g. [7, Lemma 4.8]. Combining both equalities, we obtain the desired result. □ 

As in the Gaussian case [7], the operator A is invariant under rotations about the origin in 
the plane M 2 . It is thus natural to work in polar coordinates (r, 9), and to develop the vorticity 
uj(r cos 9, r sin 9) in Fourier series with respect to the angular variable 9. In these variables, the 
action of A can be described fairly explicitly. Let 

<Kr) = y£, g(r ) = r>0. (2.21) 

Then we have the following result : 
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Lemma 2.2 Fix n G N. If uj = a n (r)s'm(n9), then Auj = n[(f>{r)a n (r) — g(r) A n (r)] cos (n6), 
where A n is the regular solution of the differential equation 

1 n 2 

-<(r) - -A' n (r) + ^A n (r) = a n (r) , r > . (2.22) 

Similarly, if uj = —a n {r) cos(n6), then Auj = n[<f)(r)a n {r) — g(r)A n (r)] sin(n6>). 

Proof. If n = 0, namely if uj is radially symmetric, it is straightforward to verify that Auj = 0. 
Thus we assume that n > 1, and that uj = a n (r) sin(n#). If i/> denotes the stream function 
defined by — Aip = uj, we have ip = A n (r)sm(n9), where A n is the regular solution of (2.22), 
namely 

An ^ = 2n (/ (r) sa ™( s ) ds + J (~) sa n( s ) ds ^) » 
The velocity field u = K[uj] = — V^V is thus 



r > . 



u = — A n (r) cos(n0)e r — A' (r) sm(n9)eg , where e r = , eg = ^— . 
r \x\ \x\ 

Since Auj = • Vuj + u • Vio» , we conclude that 

Auj = — — 7T- na n (r) cos (nO) H — A n (r) cos (n6*) — </(^ 2 ) , 

which, in view of (2.21), is the desired result. The case where uj = —a n (r) cos(n6>) is similar. □ 

As an application of Lemma 2.2, we can characterize the kernel of the operator A. We 
already know that Auj = if uj is radially symmetric. Moreover, differentiating the identity 
• Vro t = with respect to x\ and X2, we obtain A(diw*) = A(d2W*) = 0. As in the Gaussian 
case [14], we conclude: 

Lemma 2.3 Ker(A) = Xq © {a\diw* + c^c^w* | ai,a2 G M}. 

Proof. Since the decomposition (2.19) is invariant under the action of A, it is sufficient to 
characterize the kernel in each subspace X n . The case n = is trivial, because Xq C Ker(A), 
hence we assume from now on that n > 1. If uj = a n (r) sm(n9) satisfies Auj = 0, we know from 
Lemma 2.2 that (j)a n — gA n = 0. In view of (2.22), this can be written in the equivalent form 



Kir) - l -A' n (r) + - 9 ^)Mr) = , r > . (2.23) 



r 

Now, the second assumption in (2.8) means that 



r 2 q(r) 4r 2 \q'(r)\ 

SU P^TT = SU P n ( \ < 4 • 

r>o q>(r) r>0 Q(r) 

Thus, if n > 2, the "potential" term (n 2 /r 2 — g/(p) in (2.23) is positive, and since A n (r) — > 
as r — > and r — > oo, the maximum principle implies that A n = 0, hence also a n = 0. Thus 
Ker(A) n X n = {0} if n > 2. 

In the particular case n = 1, it is easy to verify that Ai(r) = rcf){r) is the regular solution 
of (2.23). Using (2.22) we find a±(r) = rg(r), so that uj = ai(r)sin(#) = —d2W*. Similarly, 
ai(r) cos(#) = — d\w*, hence the kernel of A in X\ is spanned by the functions {diw*, c^tu*}. □ 
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Using the same arguments as in [14], one can show that the operator A is not only skew- 
symmetric, but also skew- adjoint in X. This implies that Ker(A) = Ran(A)- 1 , hence 



Ran(A) = Ker(A) ± . 

Let 

Y = {uj G X | \x\ 2 uj G X} . (2.24) 

We now show that Ker(A)- 1 n Y C Ran(A), and we establish a semi-explicit formula for the 
inverse of A on that subspace. 

Proposition 2.4 If f G X n C\Y for some n>2, there exists a unique uj G X n nD(A) such that 
Auj = f. Specifically, if f = b n (r) cos(n6), then uj = a n {r) sm(n6), where 

a n {r) = p-A n (r) + , (2.25) 

(pyr) ncpyr) 

and A n is the regular solution of the differential equation 

-<M-;^> + (S-f^KM = ^> '•>»• < 2 ' 26 ' 

Similarly, if f = b n (r)sm(n6), then uj = — a n (r) cos(n6). 

Proof. If uj = a n (r) sin(n#), then Auj = n[<p(r)a n (r) — g(r)A n (r)] cos(n#) by Lemma 2.2, where 
A n satisfies (2.22). The equation we have to solve is therefore n((j)a n — gA n ) = b n , which gives 
(2.25). Moreover, combining (2.25) and (2.22), we obtain (2.26). 

Proceeding as in [5, Lemma 3.4], we now show that (2.26) has a unique regular solution, and 
we establish a representation formula. As we observed in the proof of Lemma 2.3, the "potential" 
term (n 2 /r 2 — g/(f>) in (2.23) is positive if n > 2. Let tp + , tp- be the (unique) solutions of the 
homogeneous equation (2.23) such that 

4>-(r) ~ r n as r — > , and 4>+{r) ~ r~ n as r — > oo . (2.27) 

By the maximum principle, the functions ip- are strictly monotone and linearly independent. 
The Wronskian determinant W = — satisfies W 1 + W/r = 0, hence W(r) = Innjr 

for some k > 0, and we also have 

ip-(r) ~ Kr n as r — > oo , and ip+(r) ~ nr~ n as r — > . 

With these notations, the unique regular solution of (2.26) has the following expression : 

Mr) . , +w r d. +Mr > rmm*. , *•><>. ^ 

Jo w \ s ) n n s ) Jr W(s) ncj){s) 

If / = b n (r) cos(nO) G X n , it is straightforward to verify that the function A n defined by 
(2.28) is continuous and vanishes at the origin and at infinity. Moreover, we know from (2.21) 
that 4>{r) ~ l/(2-/rr 2 ) as r — > oo. Thus, if we assume that / G X n n Y, we see that the function 
a n defined by (2.25) satisfies J °° a n {r) 2 p{r 2 )r dr < oo. As a consequence, if uj = a n {r) sin(n#), 
we conclude that uj G X n n D(X), and Auj = f by construction. □ 

Remark 2.5 If n = 1, the conclusion of Proposition 2.4 fails because djiu* G X\ n Ker(A) for 
j = 1,2. However, if f G X\C\Y satisfies (f,djW*) = for j = 1,2, one can show that there 
exists a unique uj G X\ n -D(A) n Ker(A)- 1 swc/i that Auj = /. 
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2.3 The perturbation expansion 



Equipped with the technical results of the previous section, we now go back to equation (2.14), 
which we want to solve perturbatively for large d. This equation can be written as Aw + [uj] = 
0, where A^[w] = u ■ Vu; + Before starting the calculations, we briefly explain why we 

expect to find a unique solution, under our symmetry assumptions. 

First, if ui G X satisfies (2.11)-(2.13), then ui G Ker(A)- 1 , hence oj is uniquely determined 
by Aw. Indeed, as was already observed, (2.13) means that Pqoj = 0. Moreover, it follows from 
(2.11), (2.12) that 

2 r 

(djW*,Lu) = / Xjtodx = , j = 1,2 , 

hence oj G Ker(A)- 1 by Lemma 2.3. Next, if oj and u = K[ui] have the symmetries (2.11), it is 
straightforward to verify that the nonlinearity in (2.14) satisfies N d [oj](xi, —x 2 ) = — N d [oj](xi, x 2 ) 
for all x = {x\,x 2 ) G M. 2 . This implies that PoN d [oj] = and (fliw*, N d [oj\) = 0. Moreover, we 
have {d 2 w*, N d [u]) = by construction, because this is the relation we imposed to determine the 
angular speed Q in (2.6). Thus, we see that N d [oj] G Ker(A)- 1 , and if we can prove in addition 
that |x| 2 AT d [o;] G X, then Proposition 2.4 (and Remark 2.5) will imply that N d [oj] G Ran(A). 
We can therefore hope to find a unique oj G Ker(A) 1 - n D(A) such that Aoj + N d [oj] = 0. 

To begin our perturbative approach, we compute the remainder term (2.16) for oj = 0, 
namely 

R d (x) = R d [0]{x) = (y*(x + 2x d ) - Sl[w*](x + x d ) x ) • Vw t (i) , x G R 2 . (2.29) 
From (2.7), we know that 

I X -L r<x> 
v *( x ) = Trrioi 1 ~ Q(\ x \ 2 )) i where Q(r) = / q(s)ds . 

Z7T \X\ Z J r 

By assumption, the term Q(|x| 2 ) decays faster than any inverse power of \x\ as \x\ — > oo, hence 
we can neglect its contribution in our calculations. For any fixed x G R 2 we thus have 

v*(x + 2x d ) = ±¥^ + ±V{x,2x d ) + 0(±) , as d -+ oo , (2.30) 

where 

v(x y) = (E+Kt _ vL 

v y^i y ) i i it 119' 

\x + y\ z \y\ z 

Setting x = (r cos 8, r sin#), y = 2x d = (d, 0), and proceeding as in [5, Lemma 3.2], we find 

V ^,V) = ^^(-lf^ sinM) , V 2 (x,y) = -£(-!)»_ cos(n0) . (2.31) 

n=l n=l 

In particular, returning to (2.30) and using definition (2.6), we obtain 

£l[w*\ = - [ (v*) 2 (x + 2x d )w*(x)dx = + of-M , as d ->• oo . (2.32) 
u J ^2 Trd \d / 

Note that the term V(x,2x d ) in (2.30) gives no contribution to the angular velocity f2[«;*]. On 
the other hand, inserting (2.30), (2.32) into (2.29) and using the expansion (2.31) together with 
the relation Vu>* = —xg(\x\), where g is defined in (2.21), we find for x = (r cos 0, r sin 6) : 

= 9 ~B ^ ( " ir ^ s ' m(n9) + °(d*) ' as d ^ 00 • (2 - 33) 

n=2 
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Motivated by this result, we now construct inductively an approximate solution of (2.14) of 
the form 

i l 

= E jr^ (n) (*) > «(*) = E^r« (B) (^) . ( 2 - 34 ) 

n=2 n=2 

where each velocity profile is obtained from a/ n ) via the Biot-Savart law (1.3). The order 
£ of the approximation is in principle arbitrary, but the complexity of the calculations increases 
rapidly with £, and we shall restrict ourselves to £ = 4 for simplicity. Of course, we assume that 
the symmetry and normalization conditions (2.11)-(2.13) hold at each order of the approxima- 
tion. In particular, we have 

f w (n) (x)dx = f x 1 uj ( - n \x)dx = f x 2 u/ n) (x)dx = , (2.35) 
Ji2 Jm.2 Jm.2 

for all n G {2, ...,£}. In view of [6, Appendix B], this implies that the velocity field u^ n \x) 
decays at least as fast as |x| -3 when |x| — > oo. It follows that the term u(—x — 2xd) in (2.16) is 
0(d~^) as d — > oo, and will therefore not contribute to for n < 4. For the same reason, 

f2[w* + to] = ^ j ^ ({v*) 2 {x + 2x d ) - u 2 (-x - 2x d )j (w*(x) + u(x)j dx = + > 

as d — > oo. Indeed, the leading term was computed in (2.32), and we know that the 

contribution of u 2 (— x — 2x d ) is negligible. Moreover, using (2.30), (2.31), and (2.35), it is easy 
to verify that f (v*) 2 (x + 2 X( i)oj{ X ) dx = 0(d~ 5 ), as d — > oo. Summarizing, we have shown that 

R d [u](x) = (v*(x + 2x d ) - n[w*]{x + x d )^ ■ VK(x)+u;(x)) + 0(^) 

= R d (x) + (^V(x, 2x d ) - ^) • Vw(x) + O (1) , (2.36) 

as d — > oo. Similarly, the quadratic term u ■ Vuj in (2.14) satisfies 

u- Vuj = -^u (2) • W 2) +0(^5) , as d^oo. (2.37) 

It is now a straightforward task to determine the first vorticity profiles in the expansion 
(2.34). From (2.36), (2.37), we know that the nonlinearity N d [oj] = u ■ Vuj + R d [uj] in (2.14) 
satisfies, for x = (r cos 6*, r sin 9) £ M 2 , 

N d [u>](x) = ^ (^sin(2fl) - ^ sm (3fl)) , as d ->• 00 . (2.38) 

Thus, to ensure that Au + A^M = C(d~ 4 ), we must impose 

AwW + ^(-l) n r n sin(n0) = , for n = 2,3 . (2.39) 

27T 



By Proposition 2.4, Eq. (2.39) has a unique solution a/™) G X n Pi -D(A) of the form 

w ( n )(x) = a„(r)cos(ra<9) , u (n) (x) = --A»(r) sm(n6)e r - A' n (r) cos(n9)e e , (2.40) 



where a n (r), A n (r) are given by (2.25), (2.26) with b n (r) = (—l) n r n g(r)/(2ir). As is easily 
verified, the symmetry conditions (2.11)-(2.13), are satisfied by the profiles for n = 2,3, 

and the velocity |u^(x)| decays like |x|~ n_1 as |x| — > 00. 
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Computing the profiles oj^\u^ is more cumbersome, but also more representative of what 
happens in the general case. First of all, the quadratic term (2.37) is no longer negligible, and 
using (2.40) for n = 2 we find 

u {2) (x) ■ W 2) (x) = £i(r)sin(40) , where B^r) = ^A' 2 (r)a 2 (r) - A 2 (r)a 2 (r)^j . 

Note that, to ensure that ■ Vw< 2 ' G X, we need an assumption on the second derivative of 
the function q appearing in (2.7). For instance, in analogy with (2.9), one can impose 

r k a" ( r) 2 

sup ,\ < oo , for all k G N . (2.41) 
r>o \q{r)\ 

Next, we must compute the contribution of uj^ to the right-hand side of (2.36). Using (2.31), 
(2.40), we obtain 

j2 _L 

(— V(x,2x d ) - —) ■Vw (2) (x) = £ 2 (r)sin(4#) + C(r)sm(26) , 

where 

B 2 (r) = -^(2a 2 (r) - ra' 2 (r)) , C(r) = -a 2 (r) . 

47T 7T 

Combining these results, we find instead of (2.38) : 

N d [u](x) = ^^(-i r ^ s in(^) + ^sin(4e) + ^sin(20) + o(l) , 

n=2 

as d — y oo, where B{r) = Bi(r) + B 2 (r) + r A g(r) / (2it). Therefore, in addition to (2.39), we must 
impose 

AJ 4) + B(r) sin(40) + C(r) sin(2#) = . (2.42) 

Using again Proposition 2.4, we see that (2.42) has a unique solution a/ 4 ** G (A4 + X 2 ) n D(A) 
of the form oj^(x) = a^r) cos(4#) + a 2 (r) cos(2#), where a±(r) is given by (2.25), (2.26) with 
n = 4 and 64 (r) = -B(r), while a 2 (r) is given by the same relations with n = 2 and b 2 (r) = C(r). 
An explicit expression of the velocity profile can also be obtained, as in (2.40). 

Summarizing, we have shown : 

Proposition 2.6 Let w* be a radially symmetric vorticity profile of the form (2.7), where the 
function q satisfies (2.8), (2.9), (2.41), and let 



w(x) = w*{x) + £ i" (n) (*) > < x ) = v *(x) + E ^ uH W > ( 2 - 43 ) 



n=2 n=2 



w/iere i/ie vorticity profiles oj^> elfl -D(A) n Ker(A)- 1 satisfy (2.39), (2.42), and the velocity 
profiles are obtained by the Biot-Savart law (1.3). Then w is an asymptotic solution of 
Eqs. (2.5), (2.6) in the sense that 

(v( x ) _ v (- x _ 2x d ) - n[w]{x + x d )^ ■ Vw(x) = oQ^j , (2.44) 

in the topology of X and uniformly onM 2 , as d — > 00. 
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The asymptotic expansion (2.34) is very natural, and it is clear that it can be performed to 
any finite order i 6 N if we make appropriate assumptions on the derivatives of the profile q, as 
in (2.41). As was already mentioned, we also conjecture that there exists an exact solution of 
(2.5) for d 3> 1 which coincides with (2.43) up to corrections of order 0(d~ 5 ). It is important to 
notice that, under the symmetry and normalization conditions (2.2)— (2.4), the exact solution (if 
it exists) and the asymptotic expansion (2.43) are uniquely determined by the limiting profile 
w*. 

Remark 2.7 Rather strong assumptions on the limiting profile were made in this section 
to ensure that the asymptotic expansion (2.43) holds in the function space X defined in (2.17), 
which we believe is naturally associated to the problem. This does not restrict the scope of our 
results here, because these conditions are automatically fulfilled by the Gaussian profiles created 
by the Navier-Stokes evolution. However, within the framework of Euler's equation, it is certainly 
interesting to construct interacting vortex pairs with more general profiles, including compactly 
supported ones. If we do not insist on controlling our expansion in the space X, the calculations 
presented in this section show that the assumptions on the function q can be considerably relaxed. 
The most important point is that Eq. (2.26) should have a unique solution for n > 2, and for 
n = 1 if the right-hand side satisfies some orthogonality conditions. This is definitely the case if 
the second inequality in (2.8) holds, but that condition does not imply that q is strictly decreasing 
and can well be satisfied if q is compactly supported. 

To conclude this section, we indicate how approximate solutions of (1.13) of the form (1.15) 
can be obtained from Proposition 2.6 by a simple rescaling. Given a, d > 0, we consider the 
situation described in (1.14) with ct\ = a.2 = a and r\ = r2 = d/2. If w* is a radially symmetric 
vorticity profile satisfying the assumptions of Proposition 2.6, we define, for all sufficiently small 
e > 0, 

4 4 

We {x) = «i t (i)+j;^ w (i) , v e ( X ) = «,w+E^ w w > ( 2 - 45 ) 

n=2 n=2 

where oj^ n \u^ are as in (2.43). Then, by construction, the vorticity distribution 

Ue {x) = ^w e (—^)+^w e (——^) , (2.46) 

where Xd = (d/2, 0), is an approximate solution of Eq. (1.13) with VL = a/{-ird 2 ). More precisely, 
it follows from (2.44) that 

d t uj e + (u e - ttx 1 ) ■ Vu e = 0(e) , as e -»• , (2.47) 

where u e is obtained from u e via the Biot-Savart law (1.3). 

3 Inviscid approximation of viscous vortex pairs 

In this final section, we describe in some detail the result of [5] in the particular case of a vortex 
pair, and we interprete it using the approximate solutions of Euler's equation constructed in 
Section 2. Given a, d > 0, we set Q = ajijxd 2 ) and we denote by ojq the vorticity distribution 
(1.14) where ct\ = 02 = a and r\ = r2 = d/2. For any v > 0, we consider the (unique) solution 
u u (x,t) of the rotating viscous vorticity equation 

d t u + (u-Qx 1 ) -Vuj = vAuj , xeR 2 , t>0, (3.1) 
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with initial data cuq. Up to a rotation of angle Qt, the vorticity distribution co u (x,t) coincides 
with the solution of the nonrotating equation (1.2) with the same initial data, which is studied 
in [5]. The advantage of using a rotating frame is that the vortex centers remain fixed, instead 
of evolving according to the point vortex dynamics (1.4). As a matter of fact, Theorem 2.1 in 
[5] establishes that u> u (-,t) — ojq as v — > 0, for any t > 0. 

To obtain a more precise convergence result, we decompose the solution of (3.1) into a sum 
of viscous vortices : 

v, ,x « „f x-xi \ a u ( x-x 2 \ 

" (x ' f) = Vt ^ \~Jrt~ ) + ~t W2 vTT^T' V ' ^ 3 - 2 ) 



where x\ = — x 2 = (d/2,0). As is shown in [5], both vorticity profiles w?i(£, t), (£, t) can be 
approximated by the same function 



fappfot) = G(0 + {^)F v (0 , ?GM 2 , t>0, (3.3) 

where G is the Gaussian profile (1.8) and the first order correction F u is constructed as follows. 
Let C be the Fokker-Planck operator 

C = + ±£ • V 5 + 1 , (el 2 , 

and A be the linearized operator (2.15) with w* = G, namely 

Aw = v G -Vw + K[w] ■ VG . 

Here we use the functional setting of Section 2 in the particular case where the asymptotic profile 
w* is the Oseen vortex G. This means that q(r) = \e~ r ^ in (2.7), and the assumptions (2.8), 
(2.9), (2.41) are clearly satisfied. With these notations, the profile F v is the unique solution of 
the linear equation 

-(1 -C)F V + KF V + A = , (3.4) 

a 

where A(£) = t^i&G^), see Section 3.3]. In polar coordinates £ = (r cos 0, r sin 9), we thus 
have 

A(0 = JLrV^sin^) = ^r 2 5 (r) sin(20) , 

where 5 is defined in (2.21). In particular, if we use the angular decomposition (2.19) of the 
function space (2.17), we see that A G X 2 , and it follows that F v G A2 too. Now, setting v = 
in (3.4), we obtain the simple equation Ai<o + A = 0, which coincides with (2.39) for n = 2. 
Since A2 n Ker(A) = {0}, we conclude that Fq = oj( 2 \ It is clear from (3.4) that the actual 
profile F v is close to F if the viscosity v is small compared to the circulation a of the vortices. 
As a matter of fact, it is shown in [5, Lemma 3.5] that 

H^-i^olU < C^— . (3.5) 
v + a 

To formulate our main approximation result, we introduce a function space with a weaker 
norm than X. Given any (3 > 0, we denote by Zp the space 

Zp = lw £ L 2 (R 2 ) \\w\\p < 00} , where \\w\\l = / \w(£)\ 2 e m d£ . 

Applying Theorem 2.5 of [5] to the particular situation considered here, we obtain 



17 



Proposition 3.1 Fix T > 0, and let u v (x,t) be the solution of the rotating viscous vorticity 
equation (3.1) with initial data loq. There exist positive constants K,/3, depending only on the 
product Q,T, such that, if uj u (x,t) is decomposed as in (3.2), then the vorticity profiles w\{^,t) 
satisfy 

max \\w\{-, t) -<p P (-, t)\\fi < K [j2) , (3-6) 
for all t G (0, T], where w^ pp (^,t) is given by (3.3). 

This result can be reformulated in a slightly different way, using the approximate solutions 
of Euler's equation constructed in Section 2. Indeed, if we set e = \pui, and if we remember 
that w* = G in the present case, we see that the inviscid profile w e defined in (2.45) satisfies 

«vs(o = + +°((|) 3/2 ) • 

Thus, combining (3.5) and (3.6), we obtain: 

Corollary 3.2 Under the assumptions of Proposition 3.1, we have 

(l/f\3/2 ij 

for all t € (0, T], where w e is the inviscid profile defined by (2.45) with w* = G. 

Finally, using the continuous inclusion Zp <^-» L 1 (IR 2 ), we can formulate an approximation 
result for the original function <jj v (x,t). If we compare (2.46), (3.2) and use the fact that G(£) 
and ^(0 are even functions of ^, we arrive at : 

Corollary 3.3 Under the assumptions of Proposition 3.1, we have 

I^KMl-^MIci. < kQ 3 '\c^^- , (3.8) 

for all t G (0, T], where u) e is the approximate steady state of the rotating Euler equation defined 
in (2.46) with = G. 

Comparing this last result with Theorem 1.1, we see that replacing a linear superposition 
of Oseen vortices by a more accurate solution of Euler's equation, which takes into account the 
deformation of the vortex cores due to mutual interaction, results in a better approximation. 
If we believe that there exists an exact solution of Euler's equation which is close to oj 6 for e 
sufficiently small, then (3.8) shows that the solution cj u (x,t) of the rotating viscous vorticity 
equation (3.1) slowly travels through a family of steady states of the inviscid equation (1.13) 
indexed by the length parameter e, which evolves diffusively according to e = \put. 

Remark 3.4 The right-hand side of (3.8) suggests that both dimensionless quantities v/a and 
vt/d 2 play an important role in the evolution of viscous vortices. It is tempting to eliminate 
one of these quantities by considering, for instance, the limit v — > while e = \put is kept fixed. 
Under the assumptions of Proposition 3.1, one may conjecture that the solution u u (x,t) of (3.1) 
satisfies, if e > is sufficiently small, 

ui u (x,e 2 /u) > oo e (x) , for all x G M 2 , 

V— K) 

where u) e (x) is an exact stationary solution of (1.13) of the form (2.46). Unfortunately, the 
results of [5] do not provide any control on co u (x,t) in the limit where Qt — > 00. 
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